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Abstract 

In an idealized way, some turbulent flows can be pictured by assemblies of 
many vortices characterized by a set of particle distribution functions. Ultra- 
sound provide an useful, nonintrusive, tool to study the spatial structure of 
vorticity in flows. This is analogous to the use of elastic neutron scattering to 
determine liquid structure. We express the dispersion relation, as well as the 
scattering cross section, of sound waves propagating in a "liquid" of identical 
vortices as a function of vortex pair correlation functions. In two dimensions, 

formal analogies with ionic liquids are pointed out. 
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I. INTRODUCTION 



Neutron scattering has been an essential experimental tool for the study of liquid struc- 
ture for the past 30 years M. The fact that, in a first Born approximation, the scattering 
cross section is directly proportional to the static structure factor has been a source both of 
data and inspiration for the elaboration and verification of current theories of liquid structure 

E- 

At first sight, a turbulent fluid bears no relation to the microscopic structure of a liquid. 
However, a turbulent fluid has vortices, and these objects couple to acoustic waves in much 
the same way as neutrons, or electromagnetic waves, couple to matter. The purpose of the 
present paper is to point out a number of theoretical developments that make this analogy 
more precise, particularly with respect to ionic liquids. 

It is well known experimentally as well as numerically [|BHT0|] that coherent vorti- 
cal structures easily appear in turbulent flows. Three-dimensional turbulent flows of high 
Reynolds number exhibit many intense, long-lived and tube-like vorticity regions of fairly 



well defined thickness and length ||10|| . These filaments have a tendency to organize parallel 
one to each other in clusters. Although they generally represent a small part of the mo- 
tion of the fluid, the filaments probably contain information on the statistics of the whole 
background. In a quite different context, numerical studies of two-dimensional decaying 
turbulence also reveal the emergence of coherent vortices with particle-like character. These 
vortices dominate the long-time motion of the fluid. Their size and number density are 
well described by scaling laws, and their dynamics are similar to the Hamiltonian motion 
of few point vortices [p]-|13|. According to these observations, two-dimensional decaying 
turbulence can be approximately described with a finite number of degrees of freedom. Fol- 
lowing this approach, statistical mechanical theories of the Euler's equation in a bounded 



domain have been developed for systems of point vortices [14-16], and extended afterwards 



to continuous distributions of vorticity [17 , |18|| . These theories where able to explain the re- 



sults of two-dimensional Hamiltonian simulations [IS], showing that point vortices of same 
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circulation sign are not randomly mixed, but have a tendency to organize in a nonuniform 
way within the domain. In infinite domains, homogeneous two-dimensional turbulence ap- 
pears to be well described by many large coherent vortices that must be spatially correlated. 
Quantitative experimental or numerical studies of these correlations are few, although two 
dimensional turbulence has become in recent years the object of controlled laboratory ex- 
periments p0| , |2l| 



Ultrasound scattering provides a powerful, nonintrusive tool to study vortical structure 
in flows. The characterization of the scattered pressure of an ultrasonic plane wave at suffi- 
ciently large distances from a vorticity distribution has motivated various theoretical studies 
fl22^p7|. In a way similar to elastic neutron or X-ray scattering in liquids 0, the scatter- 
ing cross section of sound waves can be related within the first Born approximation to the 
modulus of the Fourier transform of the vorticity field. This result has been experimentally 
checked for regular laboratory flows p8| , p9| , and there is increasing interest in applying this 



method for the study of turbulence f30H32f. In a related development, time-reversal acoustic 



mirrors have been used to probe vortical flows 

In the present article, we analyze this scattering problem from a distribution function 
point of view, picturing the vorticity distribution as an assembly of many undeformable 
vortices whose structure is characterized a priori by a given set of particle distribution 
densities. We will consider not only the scattering cross section problem, but also the 
dispersion relation (or effective wave number) of an acoustic wave traveling through such a 
medium. We consider a plane wave that propagates in an infinite statistically homogeneous 
flow composed of vortices with number density n. According to multiple scattering theory 
(see e.g. refs. [[34j-|36|]), such a medium can be described on average by an effective index 



of refraction, which generally depends on the wave frequency. Moreover, the wavenumber 
that characterizes the mean wave has an imaginary part, equal to half the total scattering 
cross section per unit volume; this accounts for the attenuation of the wave amplitude due 
to the loss of coherence during the scattering processes. In order to get as much insight as 
possible while keeping algebraic complications to a minimum we take all vortices with the 



same absolute value of circulation. 

In a recent publication ]37| we derived a general dispersion relation for a sound wave 



propagating in a disordered flow of mean velocity zero, and this result was applied to a 
"gas" of statistically independent identical vortices. In the next section, we extend the 
expression of the dispersion relation for a population of identical vortices of arbitrary spatial 
correlations. In section III, we specify the study to two dimensional systems and underline 
some formal analogies with ionic liquids. The results are discussed in section IV, especially 
in connection with two-dimensional decaying turbulence. Technical details are given in the 
Appendix. 

II. "LIQUIDS" OF VORTICES 

The treatment of the interaction between a sound wave and a flow presented in ref. 



37| follows an usual scheme, already presented elsewhere [26]. The analysis is based on a 
linearization of the Euler and continuity equations for a fluid, and relies on several assump- 
tions: the velocities associated with the acoustic wave are supposed to be much smaller than 
the typical velocities of the base flow, themselves much lower than the speed of sound. In 
this case the acoustic wave is a first order effect over the base flow, that can be considered 
as incompressible to zeroth order. The frequency vq of the incoming wave is also supposed 
much higher than the frequencies associated with the flow, which is considered as frozen, 
and viscous effects can be neglected. The medium is then stationary and averages are made 
over flow configurations. In the following, we denote c the speed of sound and k = v q /cq 
the wavenumber of the acoustic wave in the flow at rest. 

If the flow is composed of N undeformable vortices, its velocity field is the sum of the 
flow fields generated by each vortex, 

JV 

u(x) = . (2.1) 

a=l 

The incident wave is scattered at the velocity inhomogeneities of the medium and when 
averages over disorder configurations are performed, the mean acoustic pressure that results 



from these scattering processes can be derived within the framework of multiple scattering 
theory [p^- |36| , |38| , |39| . The mean wave is the sum of coherent scattering paths and has the 
same direction of propagation than the incident wave, noted ko, but is described by an 
effective wavenumber k 7^ ko- For a weakly perturbed medium, i.e. for a flow of low Mach 
number M. = \u\/co (the vortex density n can be high), k can be expanded perturbatively 
around ko in multiple scattering series. Since the flow has no mean velocity, the term of 
order M. is identically zero. The first non vanishing corrections involve two diagrams of 



order M. , and the dispersion relations can be written in Fourier space as [[37 



k 2 = ki + e Q _l j ^ d Y,([h.^\mh^ b \-q)]) (2.2) 

\ / a,b 

+k 2 5 7/(^3 E([^.^H?-4)][fco.^(4-g1])(fco.g1 2 ^^g1 , 

where 7 = c p /c v is the ratio of specific heats at constant pressure and constant volume, V 
the volume occupied by the flow, d the space dimension and 

G { £{q) = lim 1 (2.3) 

is the Fourier transform of the free-space Green's function. The second term of the right- 
hand side of Eq. (|2.2|) is real and proportional to the kinetic energy of the flow per unit 
volume. The real part of the last term of the right-hand side is also of order of the kinetic 
energy per unit volume, but introduces dispersion. This last term has an imaginary part, 
and using the residue theorem, one deduces that Xm{k) = A -1 is given by 

A- 1 = f = ^^/^E ([ko-^KM ~ 4)]M 6) (£o - M))) (ko4? , (2.4) 

V / a,b 

where is the total scattering cross section per unit volume, q = q/\q\, is the solid 

angle measure, and A is then the attenuation length of the coherent wave 0]. It is assumed 
that the attenuation per wavelength is small (Ako ^> 1), a condition that is fulfilled for a 
wide range of wavelengths if the Mach number is small. (See ref. [[Il]] for a fuller discussion 
of the validity of Eqs.Q-flgg)). 
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If the vortices are all identical, the Fourier transform v^ a \q) of the flow field generated 
by any vortex (a) has the form: 

y( a \q) = J dx v(x - x a , a) e~ m , (2.5) 

where x a is the position of the vortex and a an angle describing its global orientation; v 
is a fixed velocity function that depends on the vortex shape. The terms in Eq. (|2.2j) are 
averaged over all the positions and orientations {x a ,d}. 

In ref. ||37|| , we have supposed that the cross-terms in the sums of Eqs.( |2.2j )-( [2.4D vanish. 



This situation is encountered either when the flow fields of the vortices (a) and (b) do not 
overlap (v^ a \v^ ~ 0), i.e. if the number density n = N/V is small, or if the density is high 
but the vortices are spatially uncorrelated ((v( a ' .v( b >) ~ 0, or mean-field- like assumption). 
If correlations between vortices are present, the off-diagonal terms of the Born integrals 



of Eqs. ([2.2|) - (|2.4|) do not vanish. It is thus convenient to introduce a pair correlation function 
for the rigid vortices, noted g. With a normalized angle measure, this distribution function 
is such that n 2 g(x a , x b , a, b)dx a dx b dddb represents the number of pairs per unit volume and 
per unit angle formed by two vortices in configurations {x a ,d} and {x b ,b}, respectively. 
Consequently, 

(vi(x - x a , d)vj(x - x b , b)j = — J dx a dx h dadb g(x a , x b , a, b)vi(x - x a , d)vj(x - x b , b). (2.6) 

For homogeneous systems with vanishing correlations at large distances, g(x a ,x b ,d,b) = 
g(x a — x b , a, b), and g — > 1 when \x a — x b \ — > oo. Using the fact that / dx a Vi(x — x a ) = 0, 
expression fl2.6p remains unchanged if g is replaced by the total correlation function h = g — 1, 
which has the advantage of a well-defined Fourier transform. Using the identity ( p.6|) , the 



dispersion relation ( p.2|) can be recast as: 



k 2 = kl + — — I j^-dddb [k £(q,a)][k .v(-q,b)] \n5(a-b) + n 2 h(q,a,b)\ (2.7) 
Cq J {2ir) ' 

+ ~T i TT^d dhd b [ko-v(q- k ,a))[k .v(ko - q,b)](k -q) 2 Qko (?) 

Cq J {IH) 

^nd(d — b) + n 2 h(q — k Q , a, b)\ , 



where v(q, a) is the Fourier transform of the flow field of the vortex located at the origin, 
v(q, a) = J dx v(x, a) exp(-iq.x). In the same way, we get from Eq.(|2.4|): 

{n5(a — b) + n 2 h(k q — k , a, 6)| . 



°t = / 77T\d d&dh [ko-v(k q - k ,a))[k .v(k - k q,b)]{k .q) 2 (2.8) 

C J {ZTT) 



q) When the vortices are uncorrelated, h = 0, we recover the results of [37,41 



III. TWO DIMENSIONAL CASE 

Two-dimensional turbulence has been the subject of many studies because of its possible 
applications in meteorology and oceanography, but also because it is the most accessible 
dimension for computational and theoretical approaches. In the present context, the two- 
dimensional study of the effects of vortex correlations on sound propagation is clearly simpler 
because of the reduced number of degrees of freedom involved. However, we hope that it 
can provide at least some first answers qualitatively valid in any dimension, besides the 
fact that the structure of turbulence itself deeply changes with the space dimensionality. 
In two dimensions the vorticity uj points in the perpendicular direction, along the z-axis. 
We consider a system composed of axisymmetric vortices, where each vortex produces an 
azimuthal velocity field (see ref. for an example) around a vortex core of circulation 



modulus T. Suppose that each vortex has a probability x + to have a circulation +r, and 
a probability x_ to have circulation — r {x + + X- = 1). The two-particle distributions are 

described by introducing three pair correlation functions, h ++ , h and , where the sign 

of the subscript refers to the respective signs of the vortex circulations. At this point, it 
is convenient to introduce distribution functions used in the formalism of ionic liquids |42 . 
Let us consider the local vortex density, namely a "charge" density, defined as 

p^(f)=$> A (f), (3.1) 

i 

where the sum runs over the two species of number density Pi(r) and charge Zi (here, the 
orientation ±1). The vortex, or charge, structure factor Szz associated with the density 



p z {f), defined by 



Szz{q) 



(p z (q)P z (-g)) 

N 



(3.2) 



can be rewritten as 



Szz{q) = l + n x 2 + h ++ (q) + x 2 _h — (q) - 2x + x^h + _(q) 



(3.3) 



From Eq.( p.7|) , we can deduce the expression for the index of refraction defined as Af = 
c/co = 7Ze(k /k), where c denotes the new phase velocity. Noting that v(x, —z) = —v(x, z), 
and assuming that the flow is isotropic, one gets 



l- n( J „ 4) /"" dqqv 2 {q)S Z z{q) 



IttcI Jo 



(3.4) 



2 ™ ^ 
Tie 

Co 



dq 



k .v(\q- ko\) ~ S zz (\q- k \)(k .q) 2 g^(q). 



As a consequence of the structure of Eq.fl3.4p, the asymptotic results obtained in ref. ET 



for the index of refraction in the limits of short and long wavelengths, remain valid here, 
provided that one replaces any quadratic factor v 2 by v 2 Szz- 
The total scattering cross-section, in turn, is given by: 

nkn 



On 



dn^ 



k .v(\k q - k \) S Z z(\k q ~ k \)(k .q) . 



(3.5) 



Using the property q.v(q) = for incompressible fluids (a correct assumption for low Mach 
number base flows), as illustrated by the geometrical construction of Figure 1, the above 
expression can be rewritten 



Or, 



d0 



sin 2 6 cos 2 8 



2ttcq J-tt 2(1 — cos ( 



(W2-2cos#) S ZZ (k V2 - 2 COS0) . 



(3.6) 



With the same techniques presented in ref. [41], it is easy to show that <7y ~ k^ when 
k — >■ oo. If Szz(q) behaves like q 13 at small q, one gets, from Eq. fl3.6|) , 



a* T ~ /cq +/3 , k 0, 



(3.7) 



where we have used the property v(q — > 0) ~ q for an axisymmetric bounded flow. 



IV. DISCUSSION 



Relations ( |3.4|) and ( |3.5|) show that the coherent propagation of an acoustic wave through 
a population of identical vortices closely depends on their spatial structure. In two dimen- 
sions, the index of refraction and the attenuation length of the wave involve a circulation (or 
" charge" , by reference to ionic liquids) structure factor. Hence, if the vortices have a simple 
shape and generate a flow field v with well known characteristics, information on this distri- 
bution function can be deduced from the study of the acoustic properties mentioned above. 
Tractable analytical expressions have only been obtained in two dimensions. However, we 
think that they provide a qualitative understanding of the probe of vortex correlations by 
ultrasound techniques, even in higher space dimensions. 

For its similarity with two-dimensional homogeneous turbulence [j43| |, let us consider 
further the case of a neutral flow, where x + = x_ = 1/2. There are mainly two ways of 
considering the vortex structure of decaying two-dimensional turbulence, since there are two 
very different length scales in the problem. Depending on the scale of interest, one can either 
picture the flow as composed of many small, nearly point-like vortices that may form large 
structures, or consider the flow as formed by these few large coherent vortices only. One can 
assume that the size of the former vortices is of order of a dissipation scale, while the size 
of the latter is of order of the integral scale, at which energy is initially injected. 

An ideal gas of vortices is characterized by Szz(q) = 1; however, Eq. flOl) shows that it 
is equivalent, with respect to the dispersion relation, to identical ordering between identical 

and opposite vortices (h ++ + h = 2h + J). We logically check that a tendency to observe 

vortex pairs of identical circulation {h ++ + h > 2h + J) would result in an increase of 

the kinetic energy of the flow, and hence the phase change. The limit of Szz{q) at small 
wavenumbers is related to that of the total scattering cross section through relation ( |3.7| ) 
if the flow produced by one vortex is bounded in space. Szz(O) indeed represents the local 
fluctuations of the number of charges, and is of particular interest. A local electroneutrality 
assumption (Szz(fy = 0) means that the "charge" (or circulation) of a given vortex is exactly 
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canceled by the total "charge" of the vortices that surround it: at length scales of the order 
of a few vortex radii, the fluid has no net rotation. 

Let us consider the flow at small scales, i. e. made of nearly point-like vortices. Although 
not much studied in the literature, vortex pair correlations have motivated theoretical works. 
An exact analytical expression for the vorticity structure factor has been derived for a class of 



two-dimensional stationary solutions of the Navier-Stokes equation ||44|| . It is characterized 
by a Debye-Hiickel-like pair distribution, Szz(q) oc q 2 / (q 2 +k 2 ). Turbulence at high Reynolds 
number clearly exhibits quite distinct statistical properties: computer simulations rather 
show that vortices with same circulation sign have the tendency to organize in domains 
|T9| . Such systems, where electroneutrality is not locally observed, are theoretically better 
described by other approaches, like the microcanonical formulations of the statistics of two 
dimensional vortices in a bounded domain |jl4|-[T8|| . These theories predict that the one-point 
probability distribution of vortices is spatially nonuniform. To obtain this result, however, 
two-point correlation functions are usually neglected at first order (Szz — !)■ 

Two-dimensional homogeneous isotropic turbulence can be conveniently described in the 
inertial range by uniform mixture of long range correlated vortices. Notice that the two- 
dimensional energy spectrum, E(q) = q(u(q).u(—q)) / '(AirV), can be reexpressed with the 
charge structure factor as E(q) = n/(An)q v 2 (q)Szz(q)- The interpretation of the classical 
two-dimensional turbulence spectral laws E(k) ~ k~^, (/i = 3 p0| , ^5| , ^6| , /i ~ 4 [|ll],[0]) m 



terms of quasi point-like vortices (v(q) ~ q^ 1 ) leads to Szz ~ q~^ L ~ l \ This is supposed 
to be the behavior of Szz i n the inertial range. When the structure factor is steep enough 
(say, if /i is significantly larger than 1), one expects the vorticity correlations to be long 
range and positive: the spatial structure differs qualitatively from that of the ideal gas given 
by a white spectrum. The common spectral laws are compatible with preferential ordering 
between same sign vortices. 

At large scales, the flow is made of large coherent vortical structures. The pair correla- 
tions between these vortices are closely related to the properties of the energy spectrum in 
the low wavenumber limit, outside of the inertial range. For instance, if the spectrum in this 
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limit is such that E(q) ~ with // > |f48 |, the electroneutrality is local beyond the size 
of the coherent vortices. In that case, one expects that the structure of the vortex system 
should be similar to the short range structure in liquids. 
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APPENDIX A: 

For identical vortices and with the help of relation (j2.5|), the diagonal terms (a) = (b) of 
the sum 

hk .u(q)\[ko.u(-Q)\ = ^E(iko^ a Hmko^ b \-q)]) (Al) 

a,b 

give the contribution 

n J da [k .v(q, a)] [k .v(— q, a)] , (A2) 

where v(q, a) = j dx v(x, a) exp(— iq.x). The off-diagonal terms (a) ^ (b) take the form 

— —-: / dxdxdx a dx b dadb g(x a - x b , a, ^e"^^"^^^-^, ale^'^o-wf^', ^)e 4 ^' (A3) 
v v 7 

and can be recast as 

n 2 / da<i& [k .v(q, a)][k .v(— q, b)]h(q, d, b). (A4) 



Summing the terms (]AJ) and flA4), the relations Q2.2|) and (|2.4|) are transformed to Eqs.( |2.7| 
and (U). 
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FIGURES 



FIG. 1. Axisymmetric vortex in two dimensions. 



16 



